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ABSTRACT 

We present four new, high-resolution VLA radio images of the Crab nebula, taken between 2001 
February 25 and April 17. The radio images show systematic variability in the Crab's radio emission 
throughout the region near the pulsar. The principal geometry of the variable features is that of 
elliptical ripples very similar to the optical wisps. The radio wisps are seen to move systematically 
outward with projected speeds of up to 0.3 c. Comparing the new radio images to our earlier ones 
from 1998 and 2000, we show there are also more slowly moving features somewhat farther away from 
the pulsar. In particular, there is a prominent moving feature to the northwest of the pulsar which 
has a projected speed of order lO'' km s^^. Striation is seen throughout the nebula, suggesting the 
presence of wave-like disturbances propagating through the synchrotron bubble. The radio images 
were taken simultaneously with HST optical observations as part of a unique observing campaign 
to obtain simultaneous, time-resolved, high-resolution images of the Crab in different wavebands. 
Comparing the radio to the optical images, we find that the radio wisps are sometimes displaced from 
the optical ones or have no optical counterparts. We also find that some optical wisps in particular, 
the brightest optical wisps near the pulsar, do not seem to have radio counterparts. In the exterior 
of the nebula, by contrast, there is generally a good correspondence between the radio and optical 
features. 

Subject headings: ISM: Individual (Crab Nebula) — Radio Continuum: ISM — supernova remnants 



1. INTRODUCTION 

The Crab Nebula is the most easily accessible pulsar 
nebula, and thus an object of singular importance. Re- 
cent results have highlighted the remarkable nature of the 
region near the pulsar in the Crab Nebula. An image of 
this region from the Chandra X-ray observatory (Weis- 
skopf et al. 2000) showed a complex geometry consisting 
of a tilted toroidal structure with a jet along the sym- 
metry axis of the torus. A sequence of images from the 
Hubble Space Telescope {HST ) by Hester et al. (1996) 
revealed astonishing details in this region. In particular, 
there is a series of elliptical ripples in this region, which 
are usually called "wisps." These wisps have long been 
known to be variable (Lampland 1921; Scargle 1969), 
but only recent observations have shown how dynamic 
they really are. They are variable on time-scales of days, 
and move outward rapidly, with speeds up to 0.7 c (Hes- 
ter et al. 1996; Tanvir, Thomson, & Tsikarishvih 1997). 
Corresponding rapid motions are seen in the X-ray (Mori 
et al. 2002; Hester et al. 2002). In Bietenholz, Frail, & 
Hester (2001; Paper I hereafter) we showed that rapid 
motions in the wisp region were also apparent in the ra- 
dio, and that at least one wisp was moving outward with 
a projected speed of 0.24c over a period of ~ 2 months. 
Wisps are not unique to the Crab, since similar radio fea- 
tures have been seen in Vela (Bietenholz, Frail, & Hank- 
ins 1991), 3C 58 (Frail & Moffet 1993), and in B1509-58 
(Gaensler et al. 2001), although none of these have yet 
been shown to be as dynamic as those in the Crab. 

The exact nature of the wisps is not yet known, but 
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they are generally thought to be associated with a shock 
in the wind from the pulsar which powers the nebula. 
As the Crab pulsar spins down, it is losing 5 x 10^* erg 
of rotational energy per second. The bulk of this energy 
emerges from the pulsar as a highly coUimated wind of 
relativistic particles and magnetic field. The particles 
are thought to be mostly electrons and positrons, ac- 
celerated to relativistic energies, with possibly a lesser 
number of ions. In the remainder of this paper, we use 
"electrons" to refer to both electrons and positrons. This 
outflow from the pulsar is ultimately responsible for the 
synchrotron emission we observe from the body of the 
nebula. The process by which this prodigious transfer 
of energy happens, however, is not yet well understood 
(see e.g., Lyubarsky 2003; Lyutikov 2003; Spitkovsky & 
Arons, 2004; Arons 2002; Begelman 1999; Hester 1998; 
Chedia et al. 1997). 

The energy spectrum of the particles in the nebula 
can be deduced from the synchrotron emission spectrum, 
which has a power-law form extending over many decades 
of frequency, implying also power-law energy spectra of 
the electrons. The mechanism by which this particle 
spectrum is generated is still not known. In particular, 
most of the present theories of the particle acceleration 
process fail to produce a substantial number of electrons 
with relatively low energies, i.e., with Lorentz factors, 7, 
less than 10"*, which are required to produce the nebula's 
radio emission (e.g., Arons 2002; Atoyan 1999). The 
electrons with 7 < 10^ dominate by number, although 
the electrons with 7 > 10*, which produce the nebular 
optical and X-ray emission dominate the energy of the 
population. The current radio brightness of the Crab 
requires a time-averaged electron injection rate over the 
lifetime of the nebula of lO'*^"'^ s~^. 

Studying the wisps is of particular importance because 
they are the only observational window on the neutron- 
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star nebular interface. The wisps are now generally 
thought to be associated with the shock which termi- 
nates the highly coUimated outflow from the pulsar and 
randomizes it before it enters the body of the nebula. It 
is also in this region that the nebular particle spectrum 
is established. Multi- wavelength observations may there- 
fore provide observational constraints which can help an- 
swer the basic questions of how the power is transferred 
from the pulsar into the nebula and how the nebular 
particle spectrum is established. A synchronized multi- 
wavelength campaign of time- resolved observations of the 
Crab was undertaken using NRAO's VLA*, the HST and 
the Chandra X-ray satellite. 

The radio observations are of special interest because, 
as mentioned, the origin of the radio emitting electrons 
presents a particular problem. The continuity of the 
emission spectrum from the radio to the optical suggests 
that the radio-emitting electrons are produced by the 
same mechanism as the higher energy ones, and the uni- 
formity of the radio spectral index over the nebula sug- 
gests that only a single acceleration mechanism produces 
all or most of the radio-emitting electrons (Bietenholz et 
al. 1997). Understanding the acceleration process, there- 
fore, will require understanding the origin of the majority 
population of 7 < 10^ radio-emitting electrons. 

In this paper, we will discuss chiefly the radio results 
and the comparison between the radio and optical results 
of the multi-wavelength campaign mentioned above. A 
discussion of the optical and X-ray observations can be 
found in Hester et al. 2002 and Mori et al. 2002. We 
describe our observations in We show a selection of 
representative images and describe the time-evolution of 
the radio synchrotron emission in [21 In we compare 
the radio images and their time-evolution to the opti- 
cal and X-ray results, and flnally in we discuss our 
findings. 

2. OBSERVATIONS AND DATA REDUCTION 

The radio observations of the Crab Nebula were car- 
ried out with the VLA in the 5 GHz band. We observed 
approximately every two weeks between 2001 February 
and 2001 April. Table ^ shows details of our four ob- 
serving runs, all of which were done using the VLA in 
the B configuration. In all four sessions, we observed at 
frequencies of 4615 and 4885 MHz with a bandwidth of 
25 MHz, used 0556-1-238 as a phase calibrator, and cal- 
ibrated the amplitude scale using observations of both 
3C 48 and 3C 286. The total time per observing session 
was ~12 hours. We chose to observe at a spaced pair of 
frequencies within the 5 GHz band for better u-v cover- 
age, which is critical when imaging an object as extended 
as the Crab. 

Our goal was to obtain time- resolved observations, 
which would complement those in the optical with the 
HST and in the X-ray with the Chandra X-ray satellite. 
Time resolved radio observations with the VLA present 
the problem that the Crab Nebula exhibits structure at 
spatial scales from ~6' down to <1". This range of spa- 
tial scales is not well sampled by any single VLA ar- 
ray configuration, and the B array configuration we used 

The NRAO Very Large Array is a facility of the National Sci- 
ence Foundation operated under cooperative agreement by Associ- 
ated Universities, Inc. 



samples only spatial scales between ~50" and ~1'.'4 at 
5 GHz. The usual approach of recovering the large scale 
structure by additional observations using more compact 
array configurations is precluded by the fact that VLA 
configuration changes occur only once every four months. 

We devised a strategy for obtaining reliable observa- 
tions using only a single array configuration in Paper I, 
and we repeat a brief description here. Any rapid time 
evolution must occur on the smaller spatial scales, since 
a speed of c represents a proper motion of only ^^3" 
per month. Our solution to the problem of the large 
scale structure is to use maximum entropy deconvolu- 
tion (AIPS task VTESS; see Cornwell & Evans 1985, 
Cornwell 1988) by supplying a low resolution support 
(default) in the deconvolution process. The deconvolved 
image is then biased to be as close to the support as is 
allowed by the data. We chose this support to be the 
same for all four of our epochs, and we used the same 
default image that was used in Paper I (made from data 
taken in the B,C and D array configurations in 1980- 
81; Bietenholz & Kronberg, 1990, 1991). We did scale 
the support according to the known radio expansion of 
0.13% year^^ (Bietenholz et al., 1991), although the ef- 
fect is small over our observing period. The use of the 
same default for each session serves to make any dif- 
ferences between our images at different epochs be only 
those that are demanded by the data. 

We also compare the images from 2001 to our earlier 
images from Paper I, taken on 1998 August 10, October 
13, and 2000 February 11, in order to determine whether 
there are any longer timescale variations. The reduction 
of these earlier data was essentially the same as for the 
present data^, and is more fully described in Paper I. We 
re-imaged the data using the same grid spacing as was 
used for the 2001 data. 

Our process of deconvolving with a default image 
causes the structure in our radio images on the large spa- 
tial scales, i.e. those >50", to be derived predominately 
from the support image, and thus from considerably ear- 
lier multi-configuration data. This larger-scale structure 
is thus not directly comparable to that in the optical and 
X-ray observations. However, as we argued above and in 
Paper I, there is little reason to expect rapid variations 
on these scales. Our chief interest in this paper is in the 
smaller scale structure. On spatial scales smaller than 
^^20" our u-v coverage is excellent, and we believe that 
our images reliably indicate the structure of the nebula, 
and in particular, changes therein from one epoch to an- 
other. To further ensure the reliability of the images we 
present, we work for the remainder of this paper only 
with spatially high-pass filtered radio images. We con- 
servatively filtered with a Gaussian of full width at half- 
maximum (FWHM) 14", which serves to isolate those 
spatial features for which we have good u-v coverage and 
thus reliable imaging. The high-pass filtering will ensure 
that our images at each epoch contain only structure re- 

^ We note that for the 2000 February 11 epoch, wo used the 
image from 1998 October 13 rather than our usual default image. 
The reason for this was that this data set had slightly different 
u-v coverage, and the deconvolution errors appeared to be larger. 
Using the 1998 October 13 image as a default reduced the spurious 
differences from the other images. These spurious differences are 
relatively easy to distinguish from the true changes, so our results 
are not dependent on the exact choice of the default image. 
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liably determined by the data and not derived from the 
support, and will reduce any spurious differences between 
the images due to possible deconvolution errors. 

The optical observations were carried out with the HST 
Wide Field Planetary Camera 2 (WFPC2) using a rel- 
atively line- free filter centered near 5500 A (F547M). 
They consisted of 24 visits between 2000 August and 
2001 April at 11 day intervals. The radio observations 
occurred during the period between HST visits 19 to 
24. The X-ray observations were carried out using the 
Chandra Advanced CCD Imaging Spectrometer (ACIS). 
The optical and X-ray observations are more fully de- 
scribed in Mori et al. 2002 and Hester et al. 2002. 

3. RESULTS 

In Figure^] we show our latest radio image of the Crab 
Nebula, taken on 16 April 2001. The size of the restoring 
beam was 1'.'4 FWHM, which was conservatively cho- 
sen as a common size for all seven epochs. The peak 
brightness was 49 mJy beam^^, and the rms background 
level was 48 /xJy beam~^ before applying the primary 
beam correction, which causes the rms background to 
rise with the distance from the center of the image, to 
~65 /j,Jy beam~^ near the edge of the nebula. Tabled 
gives the peak brightness and rms background for all our 
images. 

Our images at the other epochs in 2001 are very similar 
to the one shown in Figure ^ and wc do not reproduce 
them here (an image from 2000 can be seen in Paper I). 
To clearly show the variation, we include an animation 
accompanying Figure^ showing the variation in the cen- 
tral region over all seven epochs of radio observations. In 
order to discuss the variable features in the nebula which 
are the focus of this paper, we proceed to make difference 
images by subtracting the image at an earlier epoch from 
that at a later one. 

3.1. Variability in the Radio Images 

We first discuss the longer term variations in the ra- 
dio images. The image in Figure |2i shows the difference 
between 2001 April 17 and 1998 August 10, showing the 
changes over an interval of ~3 years (981 days)®. Fig- 
ure 13) shows the difference between the 2001 April 10 
and 2001 February 25 images, showing the changes over 
a shorter interval of ~2 months (51 days). The most 
prominent differences over both time periods are in the 
central region of the nebula and consist principally of 
ripples with an approximately elliptical geometry. These 
ripples are very similar to what we reported in Paper I 
where we first showed the elliptical geometry of the radio 
variation over intervals of ~ 2 months and ~ 1.5 years. 

With the present data we have better time-sampling 
over both longer and shorter time-scales, and can more 
reliably determine the motions over both long and short 
time-scales. The differences over the longer time inter- 
val of ~3 years are considerably larger than those over 
the shorter interval of ^2 months, especially near the 
center of the nebula. The average rms variability of the 
central region (2' x 2' centered on the pulsar) over the 
~3 year interval is 0.4 mJy beam^^ while that over the 

^ We expanded the 1998 image by a factor of 1.0035 to account 
for the overall expansion of the radio nebula of 0.13% per year 
(Bietenholz et al. 1991) before forming the difference image. 



~2 month period is only 0.2 mJy beam^^, representing 
a fractional rms variability of 1.3% and 0.7%, respec- 
tively. The observation of larger amplitude differences 
over longer time scales is not dependent on the partic- 
ular selection of images: the differences between 2001 
March 27 and 1998 October 13 are similar in amplitude 
to those in Figure while those between 1998 October 
13 and 1998 August 10 are similar to those in Figure|2Ji. 
A significant component of the variability has a timescale 
longer than the two month observing period in 2001. In 
particular, this suggests significant variability on a time 
scale of a few years, long compared to the rapid wisp 
motions, but short compared to the age of the Nebula. 

The geometry of the variability in the central region 
consists principally of a series of elliptical arcs, remark- 
ably similar to the geometry of the elliptical wisps seen 
in the optical and X-ray (Hester et al. 1995, 1996, 2002; 
Mori et al. 2002; Weisskopf et al. 2000). In general, the 
difference features seem to consist of positive ridges lying 
farther from the pulsar than parallel negative ones, im- 
plying outward motion. The animations accompanying 
Figure n show more clearly that the motion is generally 
directly outward. 

On the difference image over the ~3 year interval 
(Fig.|2i), there are additional striations visible over most 
of the body of the nebula. The geometry is still generally 
arc-like but no longer centered on the pulsar. The pat- 
tern of positive ridges lying farther from the pulsar than 
negative ones is less readily discernible in this region, but 
the most prominent ridges still show this pattern, again 
implying generally outward motion, albeit at lower ve- 
locities and amplitudes than near the pulsar. 

Also on this '^3 year difference image, there is a promi- 
nent feature which does not share the elliptical geometry 
of the wisps. It is somewhat beyond the wisp region, to 
the northwest (NW) of the pulsar, and is in roughly the 
same orientation but on the opposite side of the jet seen 
in the X-ray and optical. We will call this feature the 
"moving arc". It extends from about 44" to 68", or 1.3 
to 2.0x10^* cm, from the pulsar^. The part nearer the 
pulsar seems to be moving transversely, while the more 
distant part is moving away from the pulsar. The proper 
motions are difficult to estimate reliably because of the 
difficulty in isolating the feature in question on the im- 
ages. From the difference images, we estimate projected 
speeds of ~7500 km s~^ transversely for the end nearer 
to the pulsar and ^14,000 km s~^ radially for the end far- 
ther from the pulsar. These speeds should be regarded 
as somewhat uncertain, since there are probably changes 
in brightness that are not associated with motion in the 
moving arc. 

3.2. Uncertainties in the Difference Images 

Are these features in the difference images real? De- 
convolution errors, caused by un-sampled parts of the 
u-v plane are often larger than the image rms background 
(Briggs, Schwab, & Sramek 1999). The deconvolution 
errors are driven by the noise, so they may be different 
from one epoch to another even in the case of identi- 
cal u-v coverage. They tend, however, to be confined to 
poorly sampled areas in the u-v plane. The result in the 

^ We adopt a distance of 2.0 kpc throughout this paper, (Trimble 
& Woltjer 1971) 
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TABLE 1 

Observing Run and Image Details 



Date (midpoint) 


Length 


Midpoint 


Image Peak Brightness 


rms background^ 




(hrs) 


(days since 1998 Aug 10) 


(mjy beam~^ ) 


(/ijy beam^-') 


2001 Apr 17 


12.1 


981 


4.83 


48 


2001 Mar 27 


12.5 


960 


4.83 


73 


2001 Mar 13 


12.5 


946 


4.81 


93 


2001 Feb 25 


12.5 


930 


4.80 


77 


. . . Additional epochs from Paper I 










2000 Feb 11 


9.4 


551 


4.77 


76 


1998 Oct 13 


10.1 


65 


4.81 


82 


1998 Aug 10 


9.6 





4.70 


100 



^'The listed values are before the VLA primary beam correction. Making this correction will cause the rms background level to rise farther 
away from the pulsar, so that the true values near the edge of the nebula are ~35% higher than those listed. 



image plane is undulations which cover much of the de- 
convolved area. Compact, well-delimited features such 
as those observed near the center of the Crab Nebula are 
not likely to be due to deconvolution errors, and neither 
are the arc-like striations seen outside the central region. 
There are some undulations on a larger scale in both dif- 
ference images, however, which likely are the result of 
deconvolution errors. 

4. COMPARISON TO THE SIMULTANEOUS HST IMAGES 

We now turn to a comparison of the radio images of 
2001 to the (almost) simultaneous HST optical ones. For 
such a comparison, we must first accurately register the 
radio and optical images, which is most easily done us- 
ing the pulsar. The pulsar is identifiable on some of the 
radio images, particularly that of 2001 March 27. The 
average radio flux density of the pulsar is ~ 0.5 mJy and 
is quite variable (Moffett & Hankins 1996), consistent 
with what is visible in our images (we note that the pul- 
sar is clearly visible when an image of the rms variation 
over the radio images is formed as we do below in Fig- 
ure |3t)- We identify the brightness peak of the pulsar in 
the radio images with that in the HST images. By us- 
ing the pulsar we are ensured of an accurate registration, 
not dependent on the absolute astrometry of either the 
radio or the HST images. The radio astrometric scale 
and rotation are such that the relative positional accu- 
racy should be a few milliarcsec. The HST astrometric 
scale and rotation are also well known, with estimated 
relative positional uncertainties over the entire image of 
<0'.'l (Raggett et al. 2002). We therefore estimate that 
the radio-optical alignment is accurate to O'.'l or better. 

We show first an almost simultaneous pair of images 
of the whole region covered by the HST CCDs, which 
is smaller than the radio images. Figure |3i shows high- 
pass filtered radio image from 2001 February 25, and 
FigurelSf) shows the optical image from 2001 February 21, 
convolved to the same effective resolution and similarly 
high-pass filtered. The animation accompanying Figure|31 
shows a side-by-side comparison of the central part of 
the radio and optical images between 2001 February and 
2001 April. We also show false color overlays of the radio 
image of 2001 April 16 with the optical images of 2001 
April 17 in Figure 0] 

In the outer part of the nebula there is generally a good 
correspondence between the radio and the optical, with 
radio counterparts existing to most optical features, with 



the obvious exception of the stars in the field. Most of 
the smaller scale structure seems to be common to the 
radio and the optical. The radio image contains some 
additional filaments, not visible in the optical, which are 
presumably in front of, or behind, the optically bright 
region, since the radio nebula is larger than the optical. 

The correlation between the radio and optical bright- 
ness is not as good near the pulsar as in the outer part 
of the nebula. To illustrate this in a statistical fashion 
we determined the correlation coefficient between the op- 
tical and the radio brightness. We sampled our images 
every 0.9" in both R.A. and decl. for this comparison, 
excluding regions near the edge of the HST CCDs and 
around the bright stars, and then determined the cor- 
relation coefficient between radio and optical brightness 
in various regions. Outside the wisp region, the correla- 
tion coefficient between the radio and optical brightness 
is 0.41 (nsampic 30, 000) while inside the wisp region it 
is 0.30 (n,ampio ~ 3000). 

We first compare the general nature of the variability 
in the wisp region in the radio with that in the optical by 
examining the rms over our images between 2001 Febru- 
ary and April. In Figure El we show images of the rms 
variation with time (a) in the radio (four images) and 
(6) in the optical (six images). The pulsar is visible in 
the radio rms image as a faint point source because, as 
mentioned, its average radio flux density is variable. The 
strongest radio variability occurs ~20" to the southwest 
(SW) of the pulsar, but radio variability is seen through- 
out the central region. The strongest optical variation is 
that of the bright wisps ^10" to the northwest (NW) of 
the pulsar, and the optical variability is concentrated also 
in the region to the NW of the pulsar although the base 
of the optical jet is also visible as a region of enhanced 
variability to the SE of the pulsar. The peak in the radio 
rms image represents a fractional variability of ~3% in 
the brightness integrated along the line of sight, while 
that in the optical represents a fractional variability of 
-7%. 

4.1. Radio Optical Comparisons near the Pulsar 

We will search for a detailed correspondence between 
the radio and the optical wisps, and seek to answer the 
important question of whether the radio wisps are at 
the same locations as the optical ones, and whether they 
move with the same speed. 
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To highlight the rapidly moving wisps, we present fur- 
ther difference images, this time between pairs of suc- 
cessive images of the closely spaced observations in 2001. 
We form both radio and optical difference images, match- 
ing the times over which the differences are formed as 
closely as possible. We show pairs of radio and optical 
difference images in Figure|51 On such difference images, 
motion of a narrow feature results in a pair of parallel 
positive and negative ridges, with the positive one on the 
side toward which the motion occurs. In order to com- 
pare the precise location of radio and optical features, 
we mark the prominent positive optical difference ridges, 
i.e., locations where the optical emission has brightened, 
in both the radio and the optical images. We indicate the 
optical difference ridges in preference to the radio ones 
because they are generally more sharply defined. 

The top panel shows the difference between the 2001 
March 13 and February 25 radio images, compared with 
the difference between the March 15 and February 21 
optical images. We mark four outwardly moving optical 
wisps to the NW of the pulsar and one to the S. The 
outermost wisp in the NW has a faint radio counterpart, 
at least within our resolution of 1'.'4. The next optical 
wisp seems distinctly displaced from a radio wisp with 
a similar geometry, with an optical positive ridge lying 
~1'.'4 outside the radio one. The remaining two optical 
wisps, which are the brightest ones, have no clear radio 
counterparts. The optical wisp to the S does seem to 
have a faint radio counterpart. There are also prominent 
radio features without obvious optical counterparts, par- 
ticularly to the SW of the pulsar. 

The middle panel shows the difference between the 
2001 March 27 and March 13 radio images, compared 
with that between the March 26 and March 15 optical 
images. The pulsar is clearly visible as a positive point 
source in the radio image but not in the optical (see § • 
As in the previous set of difference images, the bright 
optical features nearest the pulsar do not have obvious 
radio counterparts. The outer, northern wisps seem to 
exhibit some displacement between the radio and the op- 
tical, including some anti-correlation where positive ra- 
dio ridges correspond to negative optical ones. There is 
also a prominent radio ridge to the SSW of the pulsar 
which does not have an optical counterpart. 

Finally, the bottom panel shows the difference between 
the 2001 April 17 and March 27 radio images, compared 
with that between the April 17 and March 26 optical im- 
ages. The pulsar has become fainter in the radio, and 
is visible as a negative source in the radio image. The 
marked wisp to the SW of the pulsar shows a clear anti- 
correlation, being negative in the radio and positive in 
the optical, that is having become fainter between April 
17 and March 27 in the radio, but having become brighter 
over almost the same interval in the optical. As in the 
previous set of difference images, the bright optical fea- 
tures near the pulsar do not have obvious radio counter- 
parts, and there are also radio features without optical 
counterparts. 

4.2. Is there a Radio Counter-Part to the Optical Jet 

A jet is visible in both the optical and the X-ray im- 
ages, and movement along it can be seen at speeds of at 
0.3 ~ 0.4c (Hester et al. 2002; Mori et al. 2002). The 
geometry of the jet in the optical is similar to that in 



the X-ray: the base of the jet is bright, it extends ap- 
proximately 1' to the SE of the pulsar, becoming wider 
and curving slightly to the S toward its end. The base 
of the jet is visible as a region of variability ~10" to the 
SE of the pulsar in the optical rms image (Figure . 
There does not seem to be an obvious radio counterpart 
to the bright base of the jet, either in the radio rms image 
(Figure |5ji) or in the images of individual epochs (e.g., 
Fig. A highly variable feature called the "sprite" 

is seen at the base of the jet in the optical (Hester et 
al. 1998). Although there is a weak, diffuse radio fea- 
ture roughly at the location of the sprite, the positional 
alignment is not exact and the radio feature is not es- 
pecially variable, so it cannot reliably be identified with 
the sprite. 

5. DISCUSSION 

We have produced time-resolved high-resolution radio 
images of the Crab Nebula which are almost simulta- 
neous with similar optical and X-ray images. We con- 
firm the existence of the rapidly moving radio wisps in 
the central 2' of Crab Nebula, first detected as spectral 
index variations by Bietenholz & Kronberg (1992), and 
subsequently identified as mobile features in Paper I. We 
show that the rapid motions of the radio wisps are defi- 
nitely directed outward, with projected speeds of < 0.3 c. 
In addition to the wisp motion, there is complex radio 
variability over longer time scales of '^3 year, including 
features farther from the pulsar which move much more 
slowly than the wisps, but still more rapidly than the 
nebular expansion. We compared the morphology and 
motions visible in the radio to those visible in the optical 
synchrotron emission. We now discuss the relevance of 
these motions in different regions of the nebula in more 
detail. 

5.1. The Wisps and Particle Acceleration in the Crab 

Nebula 

What insight can the observed behavior of the radio 
and optical wisps give into the establishment of the neb- 
ular particle spectrum and in particular the question of 
the origin of the radio-emitting electrons? We proposed 
in Paper I that the rapidly variable radio wisps are asso- 
ciated with the optical wisps and the X-ray torus. The 
fact that the radio and optical wisps show a similar over- 
all geometry of elliptical arcs, and the fact that both are 
seen to move away from the pulsar with speeds which are 
substantial fractions of c suggests that the radio wisps are 
intimately associated with the optical ones. The elliptical 
geometry suggests that the wisps are intrinsically circu- 
lar features in or near the equatorial plane of the pulsar. 
The fractional time variability in the radio brightness is 
of the same order of magnitude as that of the optical 
brightness, suggesting that the radio-optical spectrum of 
the time- variable component, i.e., the wisps, is roughly 
comparable to that of the bulk of the nebula^. This fur- 
ther implies that the particle spectrum in the wisps is 

* The fractional variability in the radio brightness is approxi- 
mately half of that in the optical. However, the radio brightness 
along lines of sight through the center of the nebula contains con- 
siderable emission from regions in front of and behind the optically 
bright region, since the radio synchrotron nebula is larger than the 
optical one. Since these optically faint regions far from the pulsar 
are not likely to be highly variable, it seems probable that the frac- 
tional time-variability in the radio would be similar to that in the 
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also similar to that in the bulk of the nebula, at least 
over the range of energies of the particles responsible for 
the radio and the optical synchrotron emission. 

While this evidence of association between the radio 
and the optical wisps is circumstantial, we feel that it 
is compelling. The simultaneous observations, however, 
show that the correspondence between radio and optical 
wisps is not exact. While some wisps are bright in both 
radio and optical, others are bright in the radio and faint 
in the optical or vice versa. This lack of detailed corre- 
spondence argues for two different but related processes 
producing the radio and the optical wisps. By contrast, 
the correspondence between the optical and the X-ray 
wisps is good (Hester et al. 2002), suggesting that the X- 
ray wisps are produced by the same process as produces 
the optical ones. 

A long-standing problem, mentioned in the introduc- 
tion, is the provenance of the radio emitting electrons, 
i.e., those with Lorentz factors of < 10^. In most the- 
oretical pictures of the Crab, the bulk Lorentz factor of 
the pulsar wind, 7w, is expected to be ~ 10^ (e.g., Wil- 
son & Rees 1978; Kennel & Coroniti 1984a, b; Melatos 
& Melrose 1996; Arons 2002). The wisps are associated 
with a shock in the equatorial region of the pulsar wind. 
It is in, or immediately downstream from, this shock 
that the nebular particle spectrum is established. The 
problem posed by the radio-emitting electrons is that it 
is difficult to produce a substantial number of electrons 
with 7 much lower than 7^ in such a shock. Thus, if 7w 
is in fact ~ 10^, then the numerically dominant radio- 
emitting electrons with 7 < 10^ are difficult to account 
for. Furthermore, if 7w ~ 10^, then the total number of 
relativistic electrons currently injected is far too small to 
account for the number of radio-emitting electrons cur- 
rently observed in the nebula, even allowing for the fact 
that the synchrotron lifetimes of the lower-energy elec- 
trons exceed the age of the nebula. 

One possible explanation for the existence of the radio- 
emitting electrons is that they are of historical origin 
(e.g., Atoyan 1999). Since the synchrotron lifetimes of 
the 7 < 10** electrons are longer than the age of the neb- 
ula, they could have been injected at some point in the 
past. A second possibility is that there may be a second 
outflow process which is responsible for producing the 
low energy electrons — perhaps the bulk of the energy 
emerges in an equatorial wind with 7^ ~ 10^, but there 
is an additional, slower, outflow, perhaps at higher lati- 
tudes, which is the source of the radio-emitting electrons. 
In neither of these cases, however, would radio emission 
currently be expected from the wisp region. Since we do 
observe radio emission from the wisp region, neither of 
these two scenarios seem likely to apply.^ 

The association of the radio and the optical wisps 
suggests rather that there are two closely related pro- 
cesses which produce the radio and the optical wisps, 
and also the low- and the high-energy relativistic elec- 

optical, i.e., on the order of 1%, if one considers similar emitting 
volumes. 

^ We note that a population of electrons with a low-energy cutoff 
such that 7 > 7inin will produce synchrotron emission with a low 
frequency tail whose brightness is oc f^/^. However, the radio 
emission from the wisps is far too bright to be merely the low 
frequency tail for any reasonable value of 7inin, for example 7 ~ 
10^. 



trons. Two candidates for a two-part process of particle 
acceleration have been put forward. In both cases, the 
high-energy particles are produced by Fermi acceleration, 
which seems capable of producing the required particle 
spectrum. Gallant et al. (2002) hypothesized that the 
pulsar wind contains ions in addition to electrons, and 
that the resonant absorption of ion cyclotron waves ac- 
celerates lower energy electrons, as was first suggested 
by Hoshino et al. (1992). Recently, Lyubarsky (2003; see 
also Lyubarsky & Kirk 2001), showed that the low en- 
ergy part of the particle population can be produced by 
rapid, forced, reconnection of the alternating magnetic 
field stripes in the termination shock of the pulsar wind. 
In this case, the pulsar wind before the termination shock 
is Foynting dominated, but a significant fraction of the 
energy is transferred to the particles immediately after 
the shock because of reconnection, with the result that 
the average particle energy grows as the particles pass 
through the shock. This would imply a bulk Lorentz fac- 
tor in the wind of only ^ 4.4x10'^ for an injection rate 
of lO""^ s^^. Some separation between the radio and op- 
tically bright wisps could be accommodated in either of 
these scenarios, but neither makes particular predictions 
as to the relative locations of the radio and optical wisps. 

The physical nature of the wisps is also a matter of 
some debate (e.g., Hester et al 1995; Chedia et al. 1997; 
Arons 1998; Lou 1998; Begelman 1999; Komissarov & 
Lyubarsky 2003; Spitkovsky & Arons 2004). Unfortu- 
nately none of the models to our knowledge directly ad- 
dress the issue of the lower-energy electrons, or make di- 
rect predictions concerning radio wisps, largely because 
of the aforementioned problem of the origin of the lower- 
energy electrons. We do not discuss the various theories 
for the origin of the wisps in detail, but list some obser- 
vational constraints provided by the observed differences 
between the radio and optical wisps. The optical wisps 
to the NW of the pulsar are considerably brighter and 
move more rapidly than those to the SE, which is plau- 
sibly explained by Doppler beaming of material moving 
out from the pulsar with at least transrelativistic veloci- 
ties, with the bright parts of the wisps to the NW being 
those where the motion is toward us. In the case of the 
radio wisps, by contrast, those to the NW of the pulsar 
are not notably brighter than those to the SE. The ob- 
served projected velocities for the radio wisps are also 
somewhat lower than those of the optical wisps. This 
suggests that the bulk flow in the radio wisps is slower. 
The somewhat more amorphous nature of the radio wisps 
may further suggest that they are somewhat more widely 
distributed above and below the pulsar's equatorial plane 
than the optical ones. In particular, the brightest optical 
wisps near the pulsar do not seem to have radio coun- 
terparts. Perhaps, the establishment of the low energy 
particle population is slower than that of the higher en- 
ergy one, and thus the number of radio emitting electrons 
is still low in the wisps closest to the pulsar. 

5.2. The Moving Arc 

The feature we have termed the moving arc, described 
in S l3.1l above. is prominent in the '^S year difference im- 
age (Figure I2I1). It had already been seen in the radio, 
showing distinct changes in brightness between two sets 
of 1.4 GHz observations in 1982 and 1987 (Bietenholz 
& Kronberg 1992). We note that the speed of its out- 
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ward motion we estimated is higher than that estimated 
by Bietenholz & Kronberg, which would imply substan- 
tial acceleration of the moving arc. However, we mea- 
sured the speed of the moving arc from the part with the 
largest changes in brightness, while Bietenholz & Kron- 
berg measured the speed of the edge toward the pulsar, 
so a direct comparison of the speeds is probably not war- 
ranted. The moving arc is also apparent in the optical: 
it is clear in an optical difference image, between obser- 
vations in 1988 and 1994 (Hester et al. 1995), with a very 
similar morphology as is seen in the radio. It probably 
also corresponds to the "dark lane" described by Oort 
& Walraven (1956) and Scargle (1969), and possibly al- 
ready observed by Lampland (1921). Its present location 
and speed is approximately consistent with that deter- 
mined by Scargle (1969), suggesting a secularly evolving 
feature moving outward through the nebula at a pro- 
jected speed of order 10^ km s~^, about ten times that 
of the nebular expansion. 

The moving arc is approximately opposite the pulsar 
from the jet seen in the optical and X-ray. Might it be 
associated with a counter-jet? Although a counter-jet 
is not clearly seen at any wavelength, there is emission 
identified as the counter-jet channel in both the opti- 
cal and the X-ray (Hester et al. 2002). The uniqueness 
and prominence of the moving arc suggest a possible as- 
sociation with the distant end of the counter-jet. The 
alignment of the moving arc with the counter-jet is not 
exact, but since the jet is seen to be curved, it seems 
likely that the counter-jet is curved also. 

A serious problem with the association of the moving- 
arc with the counterjet is the absence of an at least sim- 
ilarly prominent radio feature associated with the jet. 
There is in fact a less prominent feature visible approx- 
imately opposite the moving arc in Figure however 
its association with the jet is entirely speculative at this 
point. Possibly the prominence of the moving arc in the 
NW over any feature in the SE is due to the more rapid 
present motion of the counter-jet, since the moving arc is 
prominent only in difference images. However, any asso- 
ciation between radio features and the optical and X-ray 
jet and its implied counter-jet must remain speculative 
until future observations show a clear connection. 

5.3. Other Motions in the Outer Nebula 

On the ~3 year difference image in FigureEh-i striations 
are visible in most of the nebula, which suggest wave-like 
motions. They are not clearly discernible enough to de- 
termine the speeds, but proper motions in the range of 
0.2 ^ 1" year~^ are suggested, implying projected speeds 
of 2 10 X lO'^ km s~^. Similar motions are reported in 
the optical in a ~6 year difference image (Hester 1995). 
The morphology suggests traveling waves, perhaps re- 
flecting off the boundary of the nebula and diffracting 
from the massive filaments of thermal gas. Since the 
wisps and the jet, where the synchrotron fluid is in- 
jected into the nebula, are seen to be unsteady on short 
time-scales, it seems reasonable to suppose that wave- 
like disturbances, excited by the wisp and jet instabili- 



ties, would propagate away from the injection region in 
the synchrotron-emitting plasma, where the sound and 
Alfven speeds are relatively high. 

6. SUMMARY 

By studying time-resolved radio observations of the 
Crab Nebula, and comparing them to simultaneous and 
similarly time-resolved optical and X-ray images we have 
found the following: 

1. We confirm the existence of radio-emitting counter- 
parts to the optical wisps. We show that the radio wisps, 
like the optical ones, have arc-like geometries and are 
moving outward at speeds of up to 0.3c. 

2. The fractional variability of the brightness at the wisp 
locations in the radio is similar to that in the optical, im- 
plying that the particle spectrum in the wisps is similar 
to that in the body of the nebula. This suggests that the 
radio-emitting electrons with 7 < lO'* are accelerated in 
the wisp region along with those with higher energies. 

3. The correspondence between the radio and the optical 
wisps is not exact. In several cases the radio brightness 
seems to be anti-correlated with that in the optical. In 
particular, the brightest optical wisps nearest the pulsar 
do not seem to have bright radio counterparts. 

4. The jet seen in the optical and X-ray does not have 
an obvious radio counterpart. 

5. The lack of detailed radio-optical correspondence sug- 
gests two different but related processes for producing 
the radio and the optical wisps, and the populations of 
relativistic electrons with 7 below and above ~ 10^. 

6. There is a prominent, curved, moving radio feature 
in the NW of the pulsar that we call the moving arc. 
It's projected speed is on the order 10'' km s~^, at least 
an order of magnitude higher than those entailed by the 
general expansion of the nebula. This is a long-lived 
feature: it was detected almost 20 years ago in the radio, 
and at least 30 years ago in the optical. 

7. There is considerable variation of the radio emission 
with time-scales of >1 year, which is longer than those 
associated with the wisps. On these time-scales, the frac- 
tional variability of the radio brightness near the pulsar 
is -1%. 

8. Striation is seen in the radio difference images span- 
ning ~3 year, suggesting the presence of wave-like mo- 
tions in the body of the nebula, albeit with lower ampli- 
tudes than the wisps and the moving arc. The geometry 
is generally arc-like but not centered on the pulsar like 
that of the wisps. 
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Fig. 1. — The Crab at 5 GHz on 2001 April 17. The FWHM size of the restoring beam was 1"4. The peak brightness is 48.5 mjy beam"^, 
and the image has been corrected for the primary beam response of the VLA antennas. The rms background was 48 ^ijy beam~^ before 
primary beam correction. The contours are drawn at 0.75, 2, 4 and 8% of the peak, and the greyscale is labelled in mJy beam~^. Maximum 
entropy deconvolution was used with a support to recover the low spatial frequency structure (see text, for details. The accompanying 
animation shows a 61" X 70" region near the center of the nebula between 1998 August 10 and 2001 April 17 (see Table [Tl . The mpeg 
time-scale is non-linear, being compressed between 1998 and 2000. The individual frames have been high-pass filtered with a 14" FWHM 
Gaussian filter to enhance the mobile features. 



10 Bietenholz et al. 




Fig. 2. — Difference images showing the changes occurring in the radio nebula, both on time-scales of a few years and about two months. 
Moving ripples are apparent throughout the central region. The greyscale is labelled in mjy beam~^ and is the same in both panels. The 
individual images were high-pass filtered with a 14" FWHM Gaussian filter before forming the differences (see text §|2j. The FWHM of 
the restoring beam is I'/A. The compact, well-delimited features near the cent er of the nebula are real, while the undulations on larger 
scales of > 14" are likely artifacts due to deconvolution errors (see text, §1^. (a): the 2001 April 17 image minus the 1998 August 10 
image, i.e., a difference image over an interval of ~3 year or 982 days, (b): the 2001 April 17 — 2001 Feb 25 difference image, i.e., over an 
interval of ~2 months or 51 days. 



Crab Wisps in Radio and Optical 11 



^ -2 a 2 4 e a -10 B fs 2a 




05 34 37 36 35 34 33 32 31 30 05 34 37 35 35 34 33 32 31 30 

RIGHT ASCENSION (JSOOO) RIGHT ASCENSION (JSOOO) 



Fig. 3. — (a) The radio image from 2001 February 25 compared with (b) the optical image from 2001 February 21, both convolved to the 
same effective resolution of 1'.'4 (FWHM; shown in lower left in each panel). Both images have been high-pass filtered with a 14" FWHM 
Gaussian filter. The radio image has been blanked to correspond to the area covered by the HST CCDs on the optical image. The radio 
greyscale is labelled in mjy beam~^, and the optical one in counts. Some edge effects due to the high-pass filtering are visible on the 
optical image. The accompanying animation shows a region of 61" X 70" near the center of both the radio (left panel) and optical (right 
panel) images between 2001 February and April (see Table IT] . The time-scale is linear. 
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Fig. 4. — An overlay of the radio image of 2001 March 27 and the optical one of 2001 March 26, convolved to the same effective resolution 
of 1'.'4. Both images have been high-pass filtered with a 14" FWIfM Gaussian filter. The radio image is in red and the optical image is in 
blue and green. Some edge effects duo to the high-pass filtering are visible on the optical image. 
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Fig. 5. — The rms variation with time over the 2001 February to 2001 April observing period. Both panels show the rms variation 
with time of images convolved to I'.'A (indicated at lower left) and high-pass filtered with a Gaussian of FWHM 14". The pulsar position 
is indicated on both panels, (a) The rms variation over the four high-pass filtered 2001 radio images, with the greyscale labelled in 
/i,Iy beain^^. (b) The rms Vciriation over the optical images of 2001 February 21, March 04, 15, 26, April 6 and 17, with the greyscale 
labelled in counts. 
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Fig. 6. — The differences between pairs of consecutive radio images on the left, with the nearest corresponding optical difference images 
on the right. Both the radio and optical images have been convolved to an effective resolution of 1'.'4 FWHM (shown in the lower left in 
each panel) and high-pass filtered with a 14" FWHM Gaussian filter. The radio greyscales are labelled in mjy beam"^. The wisps are 
more well defined in the optical, and accordingly we mark the most prominent positive difference ridges in the optical image with white 
lines. The corresponding lines are drawn in the radio image to show the location of the optical features in the radio images, switching from 
white to black lines where necessary for higher contrast. 



